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I.  INTRODUCTION 


The  long,  slender,  axisymmetric  bodies  that  are  conventionally  used  for 
spin-stabilized  projectile  designs  consist  of  an  ogive  nose  section  followed 
by  a  constant  diameter  cylinder  and  a  boattail.  The  ogive  is  similar  in  shape 
to  a  cone  but  typically  has  a  smaller  associated  wave  drag  and  has  more  volume 
for  a  payload.  Secant  ogives  that  are  normally  used  are  not  tangent  to  the 
cylinder  at  the  end  of  the  nose  section.  In  order  to  decrease  base  drag,  a 
boattail  is  added.  The  boattail  is  a  conical  frustum  that  narrows  the  base  to 
typically  75%  of  the  cylindrical  body  diameter.  This  description  character¬ 
izes  the  conventional  spin-stabilized  secant-ogive  cylinder  boattail  (SOCBT) 
projectile.  Afterbody  designs  that  statically  stabilize  nonspinning  projec¬ 
tiles  are  lifting  surfaces  of  two  types:  thin  fins  and  conical  flares.  Fins 
positioned  as  a  tail  generate  normal  force  to  statically  stabilize  the  projec¬ 
tile,  whereas  conical  flared  afterbodies  generate  high  pressure  regions  and 
move  the  center  of  pressure  rearward,  behind  the  center  of  gravity. 

The  purpose  of  this  study  is  to  evaluate  the  ability  of  a  boattai 1-flare 
(see  Figure  1)  to  improve  the  stability  of  the  SOCBT  configuration.  Computa¬ 
tional  methods  are  used  to  solve  the  parabolized  Navier-Stokes  equations  in 
the  flowfield  around  the  projectile.  Predictions  of  static  pitching  moment, 
drag  force,  and  roll  damping,  Magnus  and  pitch  damping  moments  are  made  to 
define  the  effects  of  the  boattai 1-flare  afterbody  on  the  aerodynamics  of  the 
ogi ve-cyl inder  projectile. 


II.  BACKGROUND 

The  gyroscopic  stability  of  a  spin- stabilized  projectile  can  be  assessed 
by  computing  sg,  the  gyroscopic  stability  factor.1 


V 


2  ( pd/V ) 2  IJ 
lp  ly  d&  CMa 


(1) 


where; 


pd/V 

P 


p 

d 


axial  moment  of  inertia,  kg  -  m3 
transverse  moment  of  inertia,  kg  -  m3 
dimensionless  axial  spinrate 
axial  angular  velocity,  rad/sec 
air  density,  kg  -  m3 
maximum  body  diameter,  m 


7.  Murphy,  C.H .,  "Free  Flight  Motion  of  Symmetric  Missile U.S.  Army 

Ballistic  Research  Laboratory  Report  No.  1216,  Aberdeen  Proving  Ground, 
Maryland  21005,  July  1963  (AD  442757). 


airspeed,  m/sec 

static  moment  coefficient,  per  radian 


If  0  <  Sg  <  1  the  projectile  is  gyroscopical ly  unstable.  Note  that  sg  is 

inversely  proportional  to  air  density;  therefore,  projecti les  that  are  qyro¬ 
scopical  ly  stable  for  standard  atmospheric  conditions  may  not  be  stable  under 
nonstandard  temperature  and  pressure.  Due  to  this  as  well  as  other  uncer¬ 
tainties  in  computing  sg,  many  designers  have  set  1.3  as  a  lower  limit  on 
gyroscopic  stability.' 

An  additional  requirement  for  spin-stabilized  projectiles  is  that  they  be 
dynamically  stable.  The  dynamic  stability  of  a  projectile  is  defined  by 
Murphy1  as 


2(CNa  ~  V  *  2  kx  2  CMPa 
:Na  '  2  CD  '  ky  2  (CMq  +  W 


where; 


CMq  +  CMa 


=  normal  force  coefficient  slope  per  radian 
=  total  drag  force  coefficient 
=  Magnus  moment  coefficient  slope  per  radian 
=  pitch  damping  coefficient 
=  axial  radius  of  gyration 
=  transverse  radius  of  gyration 


The  range  of  values  of  Sj  for  dynamic  stability  are  also  given  by  Murphy1 
and  are  found  to  be  a  function  of  the  reciprocal  of  the  gyroscopic  stability. 

1/sg  =  sd  (2  -  s^)  (3) 


2.  "Engineering  Design  Handbook  -  Design  for  Control  of  Projectile  Flight 
Characteristics ,  Headquarters,  US  Army  Materiel  Cormand ,  September  1966. 
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III.  COMPUTATIONAL  TECHNIQUE 


Evaluation  of  the  stability  of  secant  ogive-cylinder  projectiles  with  a 
boattail,  or  boattai 1-flare  afterbody,  requires  the  calculation  of  the  aerody¬ 
namic  forces  (normal,  drag,  and  Magnus)  and  moments  (pitching,  Magnus,  and 
roll  damping).  These  parameters  were  determined  using  three  techniques:  solu¬ 
tion  of  the  parabolized  Navier-Stokes  equations,  the  method  of  characteris¬ 
tics,  and  semi -empi ri cal  models. 

Recent  papers3,4,5,6,?  describe  the  application  of  the  thin-layer  para¬ 
bolized  Navier-Stokes  (PNS)  computational  technique  to  predict  the  flow  about 
slender  bodies  of  revolution  with  sharp  noses  at  supersonic  velocities.  The 
PNS  technique  was  used  in  this  study  to  yield  the  aerodynamic  coefficients 
necessary  for  the  evaluation  of  static,  gyroscopic  and  dynamic  stability. 

Schiff  and  Sturek4  have  employed  the  PNS  computational  technique  to  solve 
for  the  flowfield  about  a  sharp  body  of  revolution  and  have  made  comparisons 
with  experimental  data.  Computed  and  measured  pressure  data  in  both  the  axial 
and  circumferential  directions  were  consistent. 

The  PNS  solutions  do  not  compute  the  base  flow  region.  Therefore,  calcu¬ 
lation  of  the  base  drag  was  performed  using  the  method  described  by  Roache  and 
Mueller.8  The  base  flow  computation  was  started  upstream  of  the  afterbody 


3.  Schiff ,  L.B. ,  and  Steger,  J.L. ,  "Numerical  Simulation  of  Steady  Supersonic 
Viscous  Flow,"  AIAA  Journal,  Vol.  18,  No .  12,  December  1980,  pp.  1421- 
1430. 

4.  Schiff,  L.B. ,  and  Sturek,  W.B.,  "Numerical  Simulation  of  Steady  Supersonic 
Flow  over  an  Ogive-Cylinder-Boat  tail  Body,"  ARBRL-TR-02363,  US  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland  21005, 
September  1981  (AD  A106060). 

5.  Sturek,  W.B.,  and  Schiff,  L.B.,  "Computations  of  the  Magnus  Effect  for 
Slender  Bodies  in  Supersonic  Flow,"  ARBRL-TR-02384,  US  Army  Ballistic 
Research  Laboratory ,  Aberdeen  Proving  Ground,  Maryland  21005,  December 
1981  (AD  A110016). 

6.  Sturek,  W.B.,  and  Mylin,  D.C.,  " Computational  Study  of  the  Magnus  Effect 
on  Boattailed  Shell ," AIAA  Journal ,  Vol.  20,  No.  10,  October  1982,  pp. 
1462-1464. 

7.  Sturek,  V.B.,  Guidos,  B.,  and  Nietubicz,  C.J.,  "Navier-Stokes  Computation¬ 
al  Study  of  the  Influence  of  Shell  Geometry  on  the  Magnus  Effect  at  Super¬ 
sonic  Speeds,"  ARBRL-TR-02501 ,  US  Army  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  Maryland  21005,  June  1983  (AD  A130630). 

8.  Mueller,  T.J.,  and  Kayser,  L.D.,  "A  Method  of  Determining  the  Turbulent 
Base  Pressure  in  Uniform  and  Non-Uniform  Supersonic  Flows, "  ARBRL-TR- 
02374,  US  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland  21005,  October  1981  (AD  A107318). 
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using  local  velocity  profiles  from  the  PNS  computation  and  yielded  the  ratio 
of  base  to  freestream  pressure  from  a  method  of  characteristics  solution.  For 
the  cylindrical  and  boattailed  afterbodies,  the  solution  was  started  at  the 
end  of  the  cylindrical  portion  of  the  projectile,  whereas  the  solution  for  the 
boattai 1 -flare  afterbody  was  started  at  the  end  of  the  flare  with  the  charac¬ 
teristic  lines  extending  into  the  base  region. 

A  second  auxiliary  calculation  was  required  for  the  pitch  damping  coeffi¬ 
cients  used  in  determining  the  dynamic  stability  factor,  sd.  (See  Equation  2.) 

A  semi -empirical  technique  developed  by  Devan  and  Mason9  computes  the  pitch 
damping  coefficient  based  on  trends  in  wind  tunnel  and  flight  test  data  for 
this  class  of  axisymmetric  bodies. 


IV.  PROCEDURE 

To  adequately  evaluate  the  effectiveness  of  different  projectile  desigr 
a  parametric  study  was  undertaken  using  the  eight  different  shapes  shown 
Figure  1  and  described  in  Table  1.  Figure  2  is  a  schematic  of  the  projecti 
shapes  studied  grouped  according  to  geometric  similarity.  All  data  generat 
in  the  investigation  are  also  grouped  in  this  fashion.  The  projectile  shap 
had  the  same  nose  section  (ogive)  and  differed  only  in  the  length  of  t  - 
cylinder  and  afterbody  sections.  Included  in  study  were  both  the  convention^, 
boattai 1  and  unconventional  boattai 1-flare  afterbody. 


TABLE  1.  PROJECTILE  PHYSICAL  PROPERTIES.  MATERIAL  DENSITY  ASSUMED 
2.77  x  10'6  kg/mm3  (ALUMINUM) 

(mm)  (Kg)  (Kg  -  mm2) 


LI  L2  L3  Mass  lx  ly 


S0C1/2 

171.6 

0 

0 

1.74 

6.38  x  102 

S0C3 

228.8 

0 

0 

2.14 

CM 

O 

»-H 

X 

O 

• 

00 

S0CBT1 

114.4 

57.2 

0 

1.65 

5.75  x  102 

9.51  x  103 

S0CBT2 

57.2 

114.4 

0 

1.40 

4.38  x  102 

• 

PO 

00 

X 

»— * 
o 

U> 

S0CBT3 

114.4 

114.4 

0 

1.81 

6.03  x  102 

1.39  x  104 

Devan,  L.,  and  Mason,  L.A.,  " Aerodynarrri.ee  of  Tactical  Weapons  to  Mach 
Humber  8  and  Angle  of  Attack  180°:  Part  II,  Computer  Program  and  Users 
Guide,"  N.S.W.C.  TR  81-358,  September  1981 . 


TABLE  1.  (continued) 


S0CBTF1 

114.4 

28.6 

28.6 

1.69 

6.02  x  102 

1.02  x  104 

S0CBTF2 

57.2 

57.2 

57.2 

1.56 

5.11  x  102 

9.48  x  103 

S0CBTF3 

114.4 

57.2 

57.2 

1.96 

6.77  x  102 

1.86  x  104 

Because  of  these  geometric  similarities  the  PNS  code  was  advanced  over 
the  3  caliber  ogive  and  onto  the  cylindrical  section  for  1  caliber.  The 
resulting  data  were  used  to  restart  the  computation  over  different  after¬ 
bodies.  In  this  way  the  computational  Mme  was  reduced  from  the  time  that 
would  be  required  if  the  PNS  solution  was  sorted  from  the  nose  of  each  of  the 
eight  bodies.  The  computations  were  carried  out  for  an  angle  of  attack  of  2° 
and  a  nondimensional  spinrate  of  0.19  (used  for  Magnus  computations). 

Using  the  Roache-Mueller  method,  the  plane  of  data  supplied  by  the  PNS 
solution  at  the  end  of  the  body  was  advanced  into  the  base  region  to  calculate 
base  pressure.  Base  pressure  data,  in  the  form  of  a  ratio  of  base  to  free- 
stream  pressure,  was  converted  to  base  drag  coefficient  form10  using: 


2  d2  (1  -  PB/PJ 

yM2 


where 


=  base  diameter  in  calibers 
=  freestream  Mach  number 


The  semi -empirical  method  of  Devan  and  Mason  was  used  to  compute  pitch 
damping  for  the  ogive-cylinder  (S0C1,  2,  3)  and  ogive-cylinder-boattail 
(S0CBT1,  2,  3)  bodies,  the  data  base  did  not  allow  accurate  computations  of 
the  boattai 1-flare  bodies  (S0CBTF1,  2,  3).  The  pitching  moment  of  a  boattail- 
flare  afterbody  was  assumed  to  be  the  same  as  that  of  the  boattailed  body  of 
equal  length  (e.g.,  S0CBTF1  and  S0CBT1).  Flares  in  general  have  been  shown  to 


10.  McCoy,  R.L.,  " McDrag  -  A  Computer  Program  for  Estimating  the  Drag  Coeffi¬ 
cients  of  Projectiles,"  ARBRL-TR-02293,  US  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland  21005,  February  1981 

(AD  A098110). 
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Increase  a  projectile's  pitch  damping;11  predictions  of  this  coefficient  are 
therefore  conservative. 

The  computed  data  were  reduced  to  the  following  aerodynamic  coefficients. 
Total  drag  was  obtained  from  the  summation  of  viscous,  pressure,  and  base  dr-ig 
coefficients.  Normal  force  slope,  C^a  (assumed  to  be  linear),  was  obtained  by 

dividing  the  normal  force  coefficient  by  .0349  radian  (2°).  Pitching  moment 
slope,  C^,  and  Magnus  moment  slope,  C^pa,  were  obtained  by  referencing  the 

moment  coefficient  to  a  point  3.6  calibers  from  the  projectile  nose  in  all 
cases,  and  dividing  by  .0349  radian.  Roll  damping  coefficient,  was 

P 

obtained  directly  from  the  PNS  computation.  Variations  of  dimensionless  spin- 
rate,  pd/V,  were  evaluated  from  Equation  1  assuming  sg  =  1.3,  while  dynamic 
stability  was  evaluated  from  Equation  2. 


V.  RESULTS  AND  DISCUSSION 


The  aerodynamic  stability  of  the  projectile  configurations  was  studied  in 
three  ways:  static,  gyroscopic,  and  dynamic  stability.  Static  stability  is 
defined  in  terms  of  the  static  pitching  moment,  gyroscopic  stability  in  terms 
of  sg  and  dimensionless  spinrate,  and  dynamic  stability  as  defined  in  Equation 
2.  a 


Typical  wall  pressure  distributions  at  Mach  3  (angle  of  attack  2°)  are 
shown  in  Figure  3  for  one  group  of  configurations,  S0C2,  S0CBT2,  S0CBTF2.  The 
effect  of  the  boattai 1-flare  afterbody  is  immediately  apparent  from  the  drama¬ 
tic  pressure  recovery  that  occurs  over  the  flared  portion  (Figure  3c).  The 
higher  pressure  on  the  windward  side  of  the  projectile  is  due  to  the  angle  of 
attack  and  is  responsible  for  the  projectile’s  normal  force.  Higher  leeward 
pressure  on  the  boattail  detracts  from  the  total  normal  force  of  the  projec¬ 
tile  (Figure  3b).  A  large  pressure  recovery  for  the  boattai 1-flare  is  shown 
in  Figure  3c. 

The  axial  distribution  of  normal  force  for  this  group  of  configurations 
is  shown  in  Figure  4.  The  local  normal  force  rises  over  the  3-caliber  ogive 
and  1-caliber  cylinder,  decreases  sharply  over  the  2-caliber  boattail  but  is 
partially  recovered  when  one  caliber  of  the  boattail  is  flared. 

The  wall  pressure  and  normal  force  recovery  effect  of  the  boattai 1-flare 
afterbody  depends  largely  on  the  value  of  freestream  Mach  number.  Using  the 
S0CBTF1  configuration  as  an  example.  Figures  5a-c  show  the  dramatic  change  in 
wall  pressure  recovery  over  the  . 5-caliber  flare.  At  Mach  2  (Figure  5a)  wall 
pressure  recovery  is  minimal  and  the  pressure  difference  from  windward  to  lee¬ 
ward  side  is  negligible.  For  Mach  3  and  4  (Figures  5b  and  5c)  the  effect  of 
the  flare  is  significant,  producing  a  large  difference  between  windward  and 


Robinson,  M.L.,  "Boundary  Layer  Effects  in  Supersonic  Floti  Over  Cylinder- 
Flare  Bodies,"  Australian  Defense  Scientific  Service  Weapons  Research 
Establishment  Report  No.  1238,  July  1974. 
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leeward  wall  pressures.  When  these  pressure  differences  are  converted  to  nor¬ 
mal  force  distributions,  the  effect  of  Mach  number  is  again  observed  (Figures 
6a-c). 

A.  Static  and  Gyroscopic  Stability 

The  ability  of  the  boattai 1-flare  afterbody  to  augment  the  normal  force 
of  the  secant  ogive  cylinder  projectile  radically  changes  the  static  pitching 
moment.  Integration  of  the  normal  force  distribution  yields  a  total  normal 
force  located  at  the  center  of  pressure.  The  static  pitching  moment  is 
defined  by  this  force  and  the  difference  between  the  center  of  pressure  and 
center  of  gravity  locations.  Comparisons  can  be  made  between  the  SOC,  SOCBT 
and  SOCBTF  configurations  on  the  basis  of  static  pitching  moment  because  this 
moment  is  referenced  to  the  same  point  for  each  body  (3.6  calibers  from  the 
nose  as  described  in  the  Section  IV).  Figure  7a,  b,  and  c  compare  the  pitch¬ 
ing  moment  characteristics  of  the  cylinder,  boattial  and  boattai 1-flare  after¬ 
bodies  in  the  three  groups  of  configurations.  (Recall  Figure  2.)  Consistent 
with  the  axial  distribution  of  normal  force,  the  SOC  and  SOCBTF  configurations 
in  each  group  exhibit  the  lowest  pitching  moment  coefficient  (significantly 
lower  than  the  SOCBT).  Comparing  all  flared  bodies,  SOCBTF,  Figure  7d  shows 
that  the  S0C6TF1  has  the  smallest  pitching  moment  coefficient  for  all  but  Mach 
4. 


Of  the  three  flared  configurations  the  S0CBTF1  has  the  shortest  boattail- 
flare  afterbody  length;  the  pitching  moment  performance  of  this  configuration 
relative  to  configurations  with  longer  flares  (e.g. ,  S0CBTF3)  must  therefore 
be  explained.  Wall  pressure  distributions  reveal  that  the  shorter  flare  main¬ 
tains  a  large  difference  between  windward  and  leeward  pressure  over  its  entire 
length  whereas  the  longer  flares  allow  the  flow  to  expand  over  their  length 
thus  losing  normal  force.  (Compare  Figure  3c  between  5  and  6  calibers  and 
Figure  5b  between  5.5  and  6  calibers.)  The  S0CBTF1  configuration  is  unique  in 
that  it  has  the  shortest  boattai 1  and  therefore  does  not  generate  the  loss  in 
normal  force  characteristic  of  long  boattai Is . (Review  Figures  3b  and  5.)  As 
the  freestream  Mach  number  increases  to  Mach  4,  the  long  flare  of  the  S0CBTF3 
generates  increasingly  more  normal  force  and  finally  overcomes  the  loss  asso¬ 
ciated  with  its  long  boattai 1.  However,  the  S0CBTF3  is  the  greatest  normal 
force  and  the  smallest  moment  arm  (cg-cp)  and  is  therefore  subjected  to  the 
smallest  pitching  moment  at  Mach  4. 

From  computed  pitching  moment  data  and  an  assumed  s„  of  1.3,  dimension¬ 
less  spinrates  (pd/V)  can  be  determined  from  Equation  1  as  described  in  Sec¬ 
tion  IV.  Figure  8a,  b  and  c  show  that  pd/V  decreases  with  freestream  Mach 
number  for  all  configurations  with  the  lowest  pd/V  associated  with  the  SOC  and 
SOCBTF.  In  each  v.ase  the  boattail-flare  afterbody  reduced  pd/V  below  that 
required  of  the  boattail  afterbody  but  not  as  low  as  the  configuration  without 
an  afterbody  (SOC).  The  shortest  boattail-flare  afterbody  (S0CBTF1),  however, 
shows  a  pd/V  comparable  to  the  corresponding  SOC  configuration  (Figure  8c). 
In  addition  the  S0CBTF1  shows  the  lowest  required  pd/V  of  all  flared  after¬ 
bodies  (Figure  8d). 


Dynamic  stability  as  defined  in  Equation  2  requires  the  aerodynamic  coef¬ 
ficients  of  drag,  normal  force,  Magnus  moment,  and  pitch  damping.  Evaluation 
of  these  coefficients  will  precede  discussion  of  the  dynamic  stability  factor, 

sd* 


Figures  9a-d  compare  the  drag  coefficients  of  all  configurations.  As 
mentioned  in  the  proceedure  section,  drag  coefficient  is  the  algebraic  sum  of 
viscous  drag  and  pressure  drag  (from  PNS  computation)  along  with  base  drag 
(from  Roache-Muel ler).  The  boattail  afterbody  configurations  (SOCBT)  have  the 
least  drag  due  to  the  narrowing  of  the  base.  This  narrowing  causes  a  flow 
expansion  that  decreases  the  base  pressure  (Pg/P^  decreases)  but  also  geomet¬ 
rically  decreases  the  base  diameter  ( dg ) .  The  combined  effect  decreases  base 
drag.  (See  Equation  4.) 


The  boattai 1-flare  afterbody  exhibits  a  significantly  larger  drag  coeffi¬ 
cient  in  all  cases  but  for  the  S0CBTF1.  The  source  of  this  drag  penalty  is 
two- fold:  the  pressure  drag  associated  with  the  oblique  shock  formed  at  the 
boattai 1-flare  junction  and  the  higher  base  drag  associated  with  the  larger 
diameter  base.  Of  all  flared  configurations,  however,  the  S0CBTF1  shows  the 
smallest  drag  rise  (Figure  9d). 


Figures  lOa-b  compare  the  Magnus  moment  coefficient  slope  of  all  config¬ 
urations.  Consider  a  spinning  axisymmetric  body  at  angle  of  attack.  The 

component  of  freestream  velocity  perpendicular  to  the  body  (cross  flow)  is 
decelerated  on  the  side  that  opposes  spin  and  is  accelerated  on  the  side  that 
is  coincident  with  spin.  This  difference  in  velocity  across  the  body  gives 
rise  to  a  difference  in  pressure  and  a  side  force  develops,  the  Magnus  force. 
This  force  becomes  more  significant  as  the  boundary  layer  thickens  on  the  aft 
end  of  the  projectile.  This  is  especially  true  of  boattail  afterbodies  where 
the  expanded  flow  accelerates  the  boundary  layer  growth.  In  contrast  a  flare 
generates  an  adverse  pressure  gradient  and  thins  the  boundary  layer.11  In 

this  study  it  was  expected  that  the  SOCBT  conf iguration  would  have  the  highest 
Magnus  moment  with  the  SOCBTF  configuration  having  the  lowest.  This  is  sup¬ 
ported  in  Figure  10.  Comparing  all  of  the  SOCBTF  configurations  the  shortest 
flare,  S0CBTF1,  had  the  smallest  Magnus  moment  (Figure  lOd). 

The  roll  damping  coefficient  data  are  illustrated  in  Figures  lla-b.  Due 
to  the  viscous  boundary  layer  a  shearing  force  is  generated  on  the  projectile 
surface  in  the  direction  opposite  the  spinrate.  (Hence,  the  coefficient  is 
usually  negative-)  Because  the  roll  damping  coefficient  arises  from  viscous 
shear,  it  is  directly  proportional  to  the  velocity  component  tangent  to  the 
body  surface.  Using  an  aerobal 1 i Stic  coordinate  system,  the  component  of 
velocity  in  the  x-direction  (along  the  body  axis,  positive  rearward)  is  desig¬ 
nated  u,  the  component  perpendicular  to  the  body  axis  is  v,  and  the  third  com¬ 
ponent  to  form  a  right-handed  axis  system  is  w.  Roll  damping  coefficient  is 
then  proportional  to  dw/dy.  As  the  radius  of  the  spinning  body  decreases, 
surface  velocity  decreases  also  (w  =  nr).  As  Figure  11  illustrates  the  SOC 
design  has  the  largest  roll  damping  because  it  has  the  greatest  body  diameter 
from  nose  to  base.  When  a  boattail  afterbody  is  used  (SOCBT)  the  roll  damping 
decreases  slightly  as  the  body  radius  decreases.  The  boattai 1-flare  afterbody 
design  increases  roll  damping  to  a  small  degree  because  it  has  an  increased 
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body  slope  —  but  the  SOCBTF  has  less  roll  damping  than  the  SOC  with  a 
constant  diameter  shape. 

From  the  discussion  above  it  is  also  clear  that  roll  damping  is  propor¬ 
tional  to  surface  area  --  the  larger  a  body  is  the  more  viscous  shear  is  gen¬ 
erated.  Figure  11a,  that  represents  the  longest  group  of  projectiles  (the 
S0C3,  S0CBT3,  S0CBTF3  are  7  calibers  long),  also  shows  the  largest  roll  damp¬ 
ing.  Figures  lib,  and  c  represent  the  two  groups  of  projectiles  that  are  6 
calibers  in  length  and  therefore  show  less  roll  damping.  The  increase  in  roll 
damping  predicted  for  the  boattai 1-flare  afterbody  could  be  viewed  as  a  disad¬ 
vantage  to  the  design  if  it  were  not  for  the  small  magnitude  of  that  increase. 
Figure  lid  illustrates  that  the  boattai 1-flare  afterbody  design  that  leads  to 
the  smallest  increase  in  roll  damping  is  not  the  shortest  flare  (which  has 
been  the  pattern  up  to  this  point)  but  the  one  with  the  shortest  cylindrical 
section,  S0CBTF2.  Reviewing  Figure  2,  it  would  seem  that  the  S0CBTF2  has  the 
smallest  overall  amount  of  surface  area. 

With  the  aerodynamic  coefficients  of  drag  and  Magnus  moment,  the  dynamic 
stability  factor,  sd,  as  given  in  Equation  2,  can  be  evaluated.  Table  2  lists 

the  values  of  sd  grouped  by  geometric  similarity.  Dynamic  stability  was 

increased  for  all  boattai 1-flare  projectiles  when  compared  to  the  straight 
cylinder  bodies.  The  boattai 1-flare  afterbody  increased  sd  in  only  a  few 

cases  when  compared  to  the  boattai  1  bodies.  At  Mach  4,  the  boattai  1-flare 
consistently  led  to  a  more  dynamically  stable  projectile.  For  the  first 
group,  the  S0CBTF1  increased  the  stability  of  the  projectiles  for  all  of  Mach 
numbers. 


TABLE  2.  DYNAMIC  STABILITY  FACTORS  (sd)* 


S0C1 

S0CBT1 

S0CBTF1 


MACH  =  2 


.273 

.353 

.360 


MACH  =  3 


.292 

.365 

.376 


MACH  =  4 


.293 

.360 

.378 


S0C2  .273 
S0CBT2  .349 
S0CBTF2  .341 


.292 

.367 

.362 


.293 

.365 

.384 


S0C3 

.307 

S0CBT3 

.427 

S0CBTF3 

.414 

.321 

.419 

.398 


.339 

.428 

.430 


♦BOXES  INDICATE  THOSE  CONFIGURATIONS  AND  MACH  NUMBERS  AT  WHICH  THE 
BOATTAIL-FLARE  AFTERBODY  HAS  INCREASED  DYNAMIC  STABILITY. 
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The  dynamic  stability  factor  involves  many  aerodynamic  coefficients,  but 
is  defined  in  a  complicated  manner.  Platou12  has  attempted  to  determine  the 
relative  importance  of  the  individual  aerodynamic  coefficients  used  with 
dynamic  stability.  Examining  s^  in  a  parametric  manner,  Platou  concluded: 

(1)  C^a  and  Cq  have  only  a  small  influence  on  s^  except  when  CNa  -  ?C(j 
is  maximum  (i.e.,  C^a  is  maximum  and  Cq  is  minimum)  s^  is  maximum. 

(2)  Sj  is  maximum  when  the  transverse  moment  of  inertia,  Iy,  is  minimum. 

(3)  sd  is  maximum  when  Magnus  moment  slope,  CMPa,  is  minimum. 

(4)  If  pitch  damping,  CMq  +  CM*  ,  is  minimum  then  s^  is  maximum.  Also, 

if  pitch  damping  is  maximum  then  the  relative  influence  of  CMPa  on 
s^  is  increased. 

The  data  computed  by  this  study  showed  the  following  trends  for  the  flare 
shapes: 


iy 

CNa 

CD 

Sip  a 

Sq  +  CMi 

smallest 

S0CBTF2 

S0CBTF2 

S0CBTF1 

S0CBTF1 

S0CBTF1 

$ 

S0CBTF1 

S0CBTF3 

S0CBTF2 

S0CBTF2 

S0CBTF2 

largest 

S0CBTF3 

S0CBTF1 

S0CBTF3 

S0CBTF3 

S0CBTF3 

Platou' s  conclusions  can  be  used  to  explain  the  findings  illustrated  in 
Table  2.  As  mentioned,  all  projectiles  showed  an  increase  in  the  dynamic 
stability  at  Mach  4.  Examining  the  data  generated  in  this  study,  it  has  been 
shown  that: 

(1)  Magnus  moment  slope  was  minimum  at  Mach  4  for  all  projectile  designs 
studied. 

(2)  Pitch  damping  was  essentially  independent  of  Mach  number. 

(3)  At  Mach  4  the  term  CNa  -  2CD  is  a  maximum. 

Also  discussed  was  the  increase  in  dynamic  stability  generated  by  the 
S0CBTF1  for  the  first  group  of  configurations  at  all  Mach  numbers.  Examining 
the  flare  data  generated  in  this  study  it  was  shown  that: 


12.  Platou ,  A.S. ,  "The  Influence  of  the  Magnus  Moment  on  the  Dynamic  Stabil¬ 
ity  of  a  Projectile ARBRL-MR-2155,  US  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland  2100b,  January  1972 
(AD  738016). 


(1)  For  all  Mach  numbers  the  S0C8TF1  has  the  highest  C^a  and  the  lowest 
Cq,  and  therefore  has  the  greatest  value  for  the  term  CNa  -  2Cq. 

(2)  For  all  Mach  numbers  the  S0CBTF1  has  the  smallest  Magnus  moment 
slope. 

(3)  For  all  Mach  numbers  the  S0CBTF1  has  the  smallest  pitch  damping. 


VI.  CONCLUSION 

An  investigation  of  the  aerodynamics  and  stability  of  a  flared  afterbody 
for  axisymmetric  projectiles  at  supersonic  speeds  has  been  completed  using 
computational  fluid  dynamic  methods  to  predict  the  aerodynamic  coefficients  of 
interest.  The  parametric  nature  of  the  study  has  lead  to  the  evaluation  of 
the  characteristics  of  the  boattai 1-flare  relative  to  the  conventional  SOCBT 
design.  Although  all  of  the  three  SOCBTF  designs  studied  can  be  characterized 
by  similar  conclusions,  the  S0CBTF1  (having  the  shortest  boattai 1  and  flare) 
has  been  shown  to  yield: 

(1)  The  smallest  static  pitching  moment  coefficient  and  required  spin- 
rate,  therefore  the  greatest  degree  of  static  stability,  of  all 
designs  with  afterbodies  (excluding  the  SOC  designs). 

(2)  The  greatest  degree  of  dynamic  stability  enhancement  of  all  design 
groups  studied  (independent  of  Mach  number). 

(3)  The  least  Magnus  moment  of  all  designs  studied  and  the  smallest 
rise  in  total  drag  and  roll  damping  of  all  designs  studied  with 
boattai 1-flare  afterbodies. 

With  the  reduction  in  the  static  instability  of  the  spin-stabilized  pro¬ 
jectile  provided  by  the  boattai 1-flare  afterbody,  the  projectile  designer  is 
afforded  a  greater  degree  of  freedom  in  choosing  values  of  spinrate,  (or  gyro¬ 
scopic  stability  factor),  axial  and  transverse  moments  of  inertia.  Even 
though  the  boattai 1-flare  afterbody  will  allow  the  designer  to  reduce  spinrate 
and  still  achieve  gyroscopic  stability  (or  hold  spinrate  constant  and  gain 
more  gyroscopic  stability)  these  degrees  of  freedom  are  not  usually  available. 
Rather,  the  projectile  designer  is  given  a  minimum  value  of  sg  and  axial  spin- 

rate  suitable  for  the  atmospheric  operating  conditions  and  the  type  of  gun 
tube  that  will  be  used  for  firing,  and  asked  to  meet  these  parameters  with  a 
projectile  that  will  also  carry  a  given  payload.  The  designer  can  then  use 
the  boattai 1-flare  afterbody  to  extend  the  range  of  values  over  which  the 
axial  and  transverse  moments  of  inertia  may  vary.  (Axial  moment  can  be 
decreased  and  the  transverse  moment  can  be  increased.) 

Although  the  present  investigation  has  established  the  aerodynamic  and 
stability  characteristics  of  a  subcaliber  flare  (boattai 1-flare)  afterbody 
for  spin-stabilized  projectiles,  subsequent  work  may  demonstrate  further 
application  to  other  ballistic  projectile  designs. 

It  has  already  been  shown  that  the  length  (transverse  moment  of  inertia) 
of  spin-stabilized  projectiles  with  the  boattai 1-flare  afterbody  can  be 
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Figure  4.  Axial  Normal  Force  Distributions;  M  =  3.0,  a  =  2° ,  S0C2,  S0CBT2 ,  S0CBTF2 
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gure  8.  Nondimensional  Spinrate  (pd/V)  vs.  Mach  No. 
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LIST  OF  SYMBOLS 


drag  coefficient 

base  drag  coefficient  (Equation  4) 
roll  damping  coefficient 

slope  of  static  pitching  moment  coefficient 

slope  of  Magnus  moment  coefficient 

pitch  damping  coefficient 

slope  of  normal  force  coefficient 

maximum  body  diameter}m 

axial  moment  of  inertia,kg  -  m3 

transverse  moment  of  inertiajkg  -  m3 

axial  radius  of  gyration 

transverse  radius  of  gyration 

cylindrical  body  length 

boattail  afterbody  length 

flare  afterbody  length 

freestream  Mach  number 

axial  angular  velocity,  rad/sec 

ratio  of  base  pressure  to  freestream  pressure 

nondimensional  spinrate  about  projectile  axis 

gyroscopic  stability  factor  (Equation  1) 

dynamic  stability  factor  (Equation  2) 

air  speed,  m/sec 

projectile  length  in  calibers 

ratio  of  specific  heats 

air  density,  kg/m3 
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